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Abstract 

We investigate a tree-level C'(a^)-accurate action, D234c, on coarse 
lattices. For the improvement terms we use tadpole-improved coefficients, 
with the tadpole contribution measured by the mean link in Landau gauge. 

We measure the hadron spectrum for quark masses near that of 
the strange quark. We find that D234c shows much better rotational 
invariance than the Sheikholeslami-Wohlert action, and that mean- link 
tadpole improvement leads to smaller finite-lattice-spacing errors than 
plaquette tadpole improvement. We obtain accurate ratios of lattice spacings 
using a convenient "Galilean quarkonium" method. 

We explore the effects of possible 0{as) changes to the improvement 
coefficients, and find that the two leading coefficients can be independently 
tuned: hadron masses are most sensitive to the clover coefficient Cp, 
while hadron dispersion relations are most sensitive to the third derivative 
coefficient C3. Preliminary non-perturbative tuning of these coefficients 
yields values that are consistent with the expected size of perturbative 
corrections. 



1 Introduction 



Lattice QCD remains the only complete implementation of nonperturbative QCD 
and so is essential for low-energy QCD phenomenology. However, simulations 
of lattice QCD rely upon brute force Monte Carlo evaluations of the QCD path 
integral, and are very costly. In recent years it has been demonstrated that this 
cost is dramatically reduced by using coarse lattices, with lattice spacings as large 
as a = 0.4fm, together with more accurate discretizations of QCD. While highly 
corrected discretizations of gluon and heavy-quark actions are now commonplace, 
less progress has been made with the much harder problem of constructing highly 
improved light-quark actions. The best light-quark actions in widespread use 
have finite-a errors proportional to a^, which are large compared with the a^, a^a^ 
errors for improved gluon actions. The problem is compounded by the fact that 
the effective lattice spacing for light quarks is 2a rather than a, because the 
light-quark action, unlike the others, involves first-order derivatives. There are 
a number of problems, like relativistic heavy-quark physics and high-momentum 
form factors, where O(a^) improvement is crucial (in particular when used in 
conjunction with anisotropic lattices) for accurate results without more or less 
uncontrolled extrapolations over large mass and/or momentum regions. In this 
paper we take a step towards remedying this situation by presenting new results 
obtained using a highly corrected lattice action for light-quarks. 

The finite-a errors can be removed, order-by-order in a, from a lattice 
lagrangian by adding correction terms: 



In principle, the coefficients q of these correction terms can be computed using 
(weak-coupling) perturbation theory, but in lattice QCD there has been a long- 
term reluctance to rely on perturbation theory for any of the ingredients in QCD 
simulations. For most of the past twenty years this has meant that no correction 
terms were included in the action, which then has only one parameter, the bare 
quark mass; the mass is tuned nonperturbatively to give correct hadron masses. 
Recently a practical technique has been developed for nonperturbatively computing 
the coefficient of the 0{a) correction|jl]]. The (9(a)-accurate quark action, originally 
discussed in 0, has led to substantial improvements over past work, but it is still 
of limited value for lattice spacings larger than 0.1-0.2 fm. 

With only a few exceptions, it is very difficult to compute the coefficients 
for O(a^) and higher corrections nonperturbatively. Thus a perturbative 
determination is the only practical alternative that permits further improvement. 
Given the advantages of very coarse lattices, we feel it is too restrictive to abandon 
perturbation theory completely. This is particularly the case since we now know 
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that perturbation theory is generally quite reliable, provided one uses tadpole- 
improved lattice operators P, ^ . In particular, perturbation theory correctly 
predicted the relatively large renormalization of the 0{a) correction to the quark 
action several years before it was confirmed in nonperturbative studies. 

The coupling constant, as{n/a), is larger on coarser lattices, and therefore 
perturbation theory is less convergent. This makes perturbation theory less 
practical for calculating such things as the overall renormalization factors relating 
lattice currents to continuum currents. The correction terms in the quark action, 
however, are suppressed by explicit powers of the lattice spacing. Consequently 
they require less high precision, and even low-order perturbation theory may suffice 
for results accurate at the few percent level. 

In this paper we derive a tadpole-improved (9(a^)-accurate quark action, 
"D234c". We compare its predictions with the those of the standard 
Sheikholeslami-Wohlert (SW) 0(a)-accurate action, and also with the original 
Wilson (W) action. To study finite-lattice-spacing errors it is not necessary to 
take the chiral limit, so we restrict our study to quark masses near the strange 
quark's mass. Since finite-a errors tyically grow with quark mass, our results 
should improve for u and d quarks. 

The important points in our analysis are: 

• We use the mean link in Landau gauge rather than the traditional plaquette 
prescription for calculating our tadpole improvement factor uq. Our reasons 
are: (1) it gives a more rotationally invariant static potential [Q; (2) it 
has been shown in NRQCD that it leads to smaller scaling errors in the 
charmonium hyperfine splitting (3) for Wilson glue, it gives a clover 
coefficient that agrees more closely with the non-perturbatively determined 
value [^]. These studies suggest that the mean-link tadpole prescription has 
smaller quantum corrections than the plaquette prescription. Of course, once 
higher order perturbative corrections are included the two prescriptions will 
come into agreement 

• After tadpole improvement, the improvement coefficients are expected to 
have quantum corrections of order which is ~ 0.4 on our coarsest 
lattice. In Monte-Carlo simulations, we systematically study the effects of 
corrections of this size, and find that the clover coefficient Cp is the only one 
whose quantum corrections will affect hadron masses significantly, and the 
third derivative coefficient C3 is the only one that affects hadron dispersion 
relations significantly. 

• We perform various non-perturbative tests of the coefficients of the 
improvement terms. We measure the hadron dispersion relation ("speed of 
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light") to check the a^A*-^-* term; vector meson (0) scahng as a check on the 
relative weight of the clover and Wilson terms; r-dependence as an additional 
check on the relative weight of the Wilson and clover terms and the effects 
of ghost branches in the quark dispersion relation. 

• We perform a rough non-perturbative tuning of the two leading coefficients of 
the D234 action, and discuss the comparison with perturbative expectations. 

• We set our overall scale from the charmonium P — S splitting. However 
for comparisons of scaling it is ratios of lattice spacings that are important, 
and we introduce a simple method for determining these more accurately, 
and with less vulnerability to systematic errors. It consists of measuring 
P — S" in a fictitious heavy quark "Galilean quarkonium" state, i.e. using an 
NRQCD heavy quark action in which relativistic corrections are not included. 
Varying the quark mass changes the size of the state, and thereby tests for 
the presence of finite-a errors. 

We have previously studied a plaquette-tadpole-improved 0(a^)-accurate 
action on isotropic lattices 0], and plaquette-tadpole-improved 0(a^)-accurate 
actions on anisotropic lattices [|, |], and found good dispersion relations and 
scaling of mass ratios. In this paper we find that mean-link tadpole improved 
D234c has the same benefits, plus much smaller finite-a errors in hadron masses. 
This is as expected, because mean-link tadpole improvement gives a larger clover 
coefficient. 



2 D234c Quark action 

Following [^], we construct a quark action that is continuum-like (at tree level) 
through 0{a^). We start with the continuum quark action: 

S = iJcM,ij„ M, = p + m, (2.1) 

If we discretize this directly, our quark dispersion relation will contain unwanted 
doublers at the edges of the Brillouin zone. To avoid this, we perform a field 
redefinition, parameterized by r, before discretizing: 

(2.2) 

where 

= = l-lrat{p-m^). (2.3) 
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Now S = ipMip, where the transformed continuum quark operator is 

M = nMcn = p + nic- \ra{p^ -ml). (2.4) 

We use p'^ = ^fi ~ ■ discretize, allowing errors of order a'^, to obtain 

the lattice D234c quark action: 



Md234c = m,(l + ^rame) + J^ItmA^'^ - y a^Ai'^ 



(2.5) 



A^T^ is the most local centered lattice discretization of the gauge-covariant n'th 



derivative g |T0|; A^^) = A^^^A^^) = A^'^^A^^\ and A^^) = A(2)A(2). The field 
strength consists of the standard clover term fIiu \ and a relative O(a^) correction 
with coefficient C2f- 

F,Ax) ^ Fl^\x) - (Af + A(^)) fI^\x) (2.6) 

At tadpole-improved tree level, all links are divided by the Landau gauge mean 
link mq, and Cp = = C^^ = C2F = 1- We will explore the effects of deviations 
from these values. 

The terms proportional to r remove the doublers from the quark dispersion 
relations, so that for generic values of r ~ 1 this is a doubler-free tree-level 0{a?)- 
accurate quark action. The derivation can be straightforwardly generalized to 
anisotropic lattices @, [l^ . 



For r = 1 there are three fairly high ghost branches in the free quark dispersion 
relation (Figure |I|). To investigate the effect of redundant terms we will also study 
r = 2/3, for which one of the ghost branches moves down so that -E(O) ~ 1.0. 

Note that the two leading terms both violate symmetries that will be restored in 
the continuum limit, and hence can be non-perturbatively tuned. The clover term 
violates chiral symmetry, and so can be tuned by imposing PCAC The A^^) 
term is the only rotational symmetry- violating 0{a?) term, and so its coefficient 
C3 can be nonperturbatively tuned by imposing rotational invariance. 
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Eip) 




E{p) 




Figure 1: Massless dispersion relation Re[ai?(a|p|)] for D234c(r = 1) and 
D234c(r = 2/3) quarks on an isotropic lattice. Note the continuum-like behavior 
at low momentum, and the presence of three ghost branches, one of which drops 
dramatically as r changes from 1 to 2/3. 



3 Gluon action and lattice spacing determination 

We use a tree-level tadpole-improved plaquette and 2x1 rectangle glue action 
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S 



X,fl 



1 Rfiiyjx) 

12 ul 



1 ^^m(^) 

12 ul 



|ReTrl4^(yU + u — ft — u) (plaquette), 
|ReTriy(2/i + z/ - 2/i - z/) (rectangle). 



(3.1) 



where a Wilson loop W{p — a ■ ■ ■) goes one link in the p direction, one link in the 
negative a direction, etc. 

This definition of jS is different from that of Ref. [|13], where a factor of 5/3^0 



was absorbed into (3. We prefer the notation here because (3 = Q/g"^ as in the 
original Wilson action. Furthermore the coupling 



a., = 



2tt(3 



(3.2) 



is now tadpole- improved and therefore roughly equal to continuum couplings like 

The tadpole improvement factor uq is the mean of the link operator in Landau 
gauge. At both our lattice spacings we found that an S'' lattice was large enough 
for finite volume effects in uq to be of order 0.1%. To fix to Landau gauge we 
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maximize 



V ^ ReTr {u,{x) - U^{x) U,{x + fi)] . (3.3) 

We generated gluon configurations at two lattice spacings, 0.40 fm and 0.25 fm 
(see Table |I|). The lattice spacing was determined in two ways. 

Firstly, we performed NRQCD simulations of charmonium, using the 
experimental value of 458 MeV for the spin-averaged P — S splitting; the results 
are given in Table |I| column 6 (details in Table |^). Note that the errors quoted are 
statistical, and do not reflect systematic uncertainties such as quenching, finite- 
a errors, or higher-order relativistic effects neglected in our NRQCD simulation. 
This lattice spacing determination is therefore not suitable for precise comparisons 
with data from other groups. 

Secondly, a more accurate determination of ratios of lattice spacings is possible, 
since there is no need to simulate a known physical state. At each lattice spacing 
we can measure the mass of a fictitious state, whose properties are chosen for 
convenience. We chose "Galilean quarkonium," a bound state of a quark and 
antiquark in a non-relativistic world. We simulated this state using NRQCD with 
no relativistic corrections. By making the Galilean quarkonium lighter (and hence 
bigger) than charmonium we reduce the finite-a errors. In fact, we studied a range 
of quark masses down to about half the charm quark mass and found that lattice 
spacing ratios were all consistent with each other, within errors. For details see 
[Appendix A | and Table ^ Since our main goal is to compare our mean-link- 
tadpole- improved results with SCRI's plaquette-tadpole-improved results, we used 
Galilean quarkonium to calculate the ratios of our lattice spacings to the SCRI 
P = 7.4: lattice spacing. The results are in Table |^. 

For convenience we want to give our results an absolute energy scale, so we 
take the SCRI /? = 7.4 lattice spacing to be = 812 MeV, which corresponds 
to ^/o^ = 468 MeV for their string tension. This gives the final column of Table |l|, 
which is consistent with our charmonium measurements. Note that the error bars 
reflect the uncertainty in the ratio to SCRI's lattice spacings, which is the relevant 
quantity for scaling comparisons. It does not reflect the overall uncertainty in the 
scale, which was introduced purely for convenience. 

For hadron spectrum measurements we used lattices of the same physical 
size (2 fm) at both lattice spacings. We also performed a set of measurements 
investigating the volume dependence of hadron masses (see appendix. Table H). 
We see that the 1.6 fm lattice agrees with the 1.75 fm and 2 fm lattices within 
statistical errors. 
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P as 




Up 


approx a 


charmonium a ^ 


rel. to SCRI 


1.157 0.413 


0.738 


0.8196 


0.40 fm 


495(4) MeV 


497(3) MeV 


1.719 0.278 


0.797 


0.8576 


0.25 fm 


790(10) MeV 


785(6) MeV 


7.4 (SCRI) 






0.24 fm 


840(20) MeV 


812(def) MeV 



Table 1: Parameters of glue for our two lattice spacings, and one of SCRI's. Mean 
Landau-gauge link (uq) was used in our tadpole improvement (TI); plaquette TI 
factor Up = (Plaq)^/^ is given for comparison. The relative to SCRI is obtained 
by the Galilean quarkonium method. 



4 Results 



We subjected the D234c action to a series of tests to determine its viability at large 
lattice spacings. We examined the scaling of the vector meson mass and of baryon 
masses, and we measured hadronic dispersion relations. We also measured the 
sensitivity of these physical quantities to changes in the coefficients in the action. 



4.1 Hadron masses 



In Figure ^ we show how the vector meson mass varies with lattice spacing when 
the ratio of the pseudoscalar to vector meson masses is P/V = 0.7. (Full data is in 
[Appendix B | , along with data for P/V = 0.76, corresponding to a slightly larger 



quark mass; see tables |10], g). We present data obtained using the tree-level D234c, 
with Landau-link tadpole improvement, at lattice spacings a of 0.25 fm and 0.4 fm. 
These values are compared with results from SCRI obtained using the Wilson and 
SW actions ^M, which can be extrapolated to give a = results, as indicated. 



We also examined baryon masses; ratios of these to the vector mass are shown 
in Fig. H Our measured values at a = 0.25 fm are within 2a of the quadratic 
extrapolation of the SCRI values to the continuum. 



a 


mean-link D234c 


mean-link SW 


plaquette SW 


plaquette Wilson 


0.4 


0.969(9) 


0.933(10) 






0.25 


1.027(10) 


1.041(15) 













1.035(5) 


1.034(37) 



Table 2: Phi (V) masses in GeV at P/V = 0.70, from tables |I| and |TU|, using a 
quadratic fit to SCRI's SW data to extrapolate it to a = 0. For the D234c action 
we see a 2(1)% scaling error at a = 0.25 fm, and 7(1)% scaling error at a = 0.4 fm. 

In table ^ we give the vector mass at our two lattice spacings, along with a naive 
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1.2 
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1 1 

D234c, mean-link Tl 
SW, mean-link Tl iG-\ 
SW, SCRI, plaq Tl \^ 
Wil, SCRI, plaq Tl hM 


1.1 










1 


1 


'x. . 




i 


0.9 






"X. 




0.8 




"i 






0.7 




i- 


1 


1 



0.05 0.1 0.15 0.2 

(fm^) 



Figure 2: Mass of meson as a function of lattice spacing, P/V = 0.70. Mean 
link Tl (• D234, o SW) clearly scales better than plaquette Tl (x SW, + Wilson, 
data kindly supplied by SCRI) 



continuum extrapolation of SCRI's data. (For details of the scale determination 
see [Appendix A j) . The errors include the relative error between our mass scale 
and SCRI's, but not the uncertainty in the absolute scale, which may be affected 
by quenching effects. 

Using a quadratic fit to SCRI's SW data to extrapolate it to a = 0, the 
D234c action shows a 2(1)% finite-a error at a = 0.25 fm, and 7(1)% finite-a 
error at a = 0.4 fm. (The quadratic fit may be too naive: the true continuum 
value could differ by a few percent, however this will not substantially affect 
our conclusions below.) These finite-a errors are due to radiative corrections 
to the tree-level coupling constants, and higher-order interactions not included 
in our action. We measured the sensitivity of the hadron masses to radiative 
corrections by varying each of the tree-level coupling constants. The fractional 
change caused by multiplying each coupling constant in turn by 1 -|- a<j is shown 
in Table |^. The only coupling for which radiative corrections are important is 
the clover coupling, Cp A perturbative analysis of Cp through 0{as) will 

soon be completed |jl6[. Alternatively, the 0{as) coefficient could be determined 
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1 1 

mean-link D234c h«— i 
plaq SW he— 1 


1.8 










D/V 


1.7 










I-*-- 

_ . .m 


1.6 










. . ■ 3E- ■ ■ ■ ' 


1.5 








...X ■■■■■■■■ $ 


1.4 








1 


N/V 

1 1 



0.05 0.1 0.15 0.2 

(fm^) 



Figure 3: D/V and N/V for D234c and plaq TI SW (SW data kindly supplied 
by SCRI.) 

by making the vector meson mass at our smallest lattice spacing agree with the 
continuum. For example, using our continuum extrapolation of SCRI's data we 
find that taking 

C^ = l + 0.2a„ (4.1) 

where is the bare coupling ( |3.2| ), reduces the V mass error at 0.25 fm from 2(1)% 
to 0(1)%, and at 0.4 fm from 7(1)% to 0(1)%. This suggests that perturbative 
corrections to Cp are relatively small after tadpole improvement. The necessary 
corrections to Cp are still perturbative if the true continuum value differs by a few 
percent from the naive extrapolation of the SCRI data. 

We also include SW results in Figure ^ using the Landau-link tadpole- improved 
tree-level value for the clover coefficient, Cp- Comparing these with SCRI's SW 
results, for which the plaquette was used to determine uq, we see that the Landau- 
link results show much smaller finite-a errors. This result was anticipated based 
on work at smaller lattice spacings using SW quarks and the (unimproved) Wilson 
action for the gluons . 
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coeff 




C p 




^4 


C2F 


a (fm) 


0.25 


0.4 


0.4 


0.4 


0.4 


P 


11(1)% 


35(1)% 


3.8(6)% 


-1.0(2)% 


-2.9(9)% 


V 


10(1)% 


35(1)% 


3.7(4)% 


-1.0(2)% 


-2.9(5)% 


N 


6(2)% 


26(2)% 


3.4(7)% 


-0.9(2)% 


-1.7(7)% 


D 


8(1)% 


25(2)% 


3.8(9)% 


-1.1(3)% 


-2.3(4)% 



Table 3: Percentage change in hadron masses when individual coefficients in the 
mean-link tree-level TI D234c quark action (Eq. |2.5|) are multiplied by 1 + a^. All 
are at P/V ^ 0.7, except the C^, a = 0.4 is at P/V = 0.76. 

4.2 Hadron dispersion relations 

The finite-a errors in the (p mass appear to be almost as small for the SW action 
as they are for D234c, but this is deceptive. To see why, we consider the quantity 

P = ^ ' 4.2 

P 

for different hadrons and three momenta p, where E{p) is the hadron's energy. In 
the continuum limit, c^(p) = 1 for all p. This quantity is particularly sensistive 
to the C3 term in the D234c action since this term is not rotationally invariant; 
it cancels the leading (non-rotationally invariant) error in the SW action. Our 
results for c^, for both pseudoscalar and vector mesons, are shown in Figure ^. 
At 0.4 fm, D234c is dramatically superior: it deviates from = 1 by only 3-5% 
at zero momentum, and by less than 10% even at momenta of order 1.5/a, while 
SW gives results that deviate by 40-60% or more for all momenta, including zero. 
As expected both formalisms improve at 0.25 fm, although D234c is still clearly 
superior. 

The results have practical implications. For example, there are two different 
definitions of a hadron's mass in lattice simulations. One is the "static mass," 

m,t.tic = E{p = 0), (4.3) 

and the other the "kinetic mass," 

mkinetic = lim = E{0)/c'{0). (4.4) 

P^o dE/dp 

In D234c these two definitions agree to within 3-5% at 0.4 fm. In SW they differ 
by 40-60%, making it impossible to say what the "true" mass is for this formalism 
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C2(p) 




0.5 



1.0 



D234c 

I I I I I I I I I I I I I I I I 

SW 



1.2 h 
1.0 
0.8 
0.6 - 
0.4 

0.2 h a = 0.25 fm 



0.5 1.0 

p (GeV) 

Figure 4: Speed of light squared using the D234c and SW actions, for the 
pseudoscalar (•) and vector (o) mesons, at P/V = 0.76, on a 5^ x 8 x 18 lattice 
at a = 0.4 fm (see Table |1^), and on a 8^ x 24 lattice at a = 0.25 fm (see Table 
131) . [D234c V points are offset for clarity.] D234c shows much better rotational 
invariance; both discretizations improve at the smaller lattice spacing, where the 
pseudoscalar and vector results for D234c are indistinguishable with our statistics. 
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at 0.4 fm. The kinetic and static masses in either formahsm must be equal for 
zero-mass mesons, because of the axis-interchange symmetry of the actions. The 
deviations seen here are because the strange quark is relatively massive at our 
lattice spacings: the mass, for example, is 2.1/a with D234c at 0.4 fm. These 
results illustrate that D234c is far more accurate than SW for hadrons with large 
masses in lattice units. More generally, D234c is far more accurate for hadrons 
with large energies and/or large momenta. 

The C3 term is the only a? correction that breaks rotational invariance. Thus 
we can use to tune C3 nonperturbatively. In our 0.4 fm simulation, we tuned C3 
to make the dispersion relation for the lightest meson, the pseudoscalar, perfect at 
low momentum: when 



Ca = 1.2 (4.5) 
^ l + 0.5a, (4.6) 

we obtain the results shown in the top part of Figure ^ The for the pseudoscalar 
is now within ±2% of c^ = l at least out to momenta of order 1.5/a. 

It is noticeable in figure || that tuning the pseudoscalar has worsened the 
dispersion relation for the vector meson (c^(p = 0) = 1.21(5) from a quartic fit). 
However, the vector is 25% heavier than the pseudoscalar and so should have 
substantially larger (amy)" errors. This is confirmed by a rerun at lower quark 
mass, P/V = 0.6, where the vector is lighter (bottom part of Figure There 
we see that chosing C3 = 1.2 works for both pseudoscalar and vector: the vector 
now has c^(p = 0) = 1.04(4). In table |lT| we see that, to within statistical errors, 
the pseudoscalar meson dispersion relation is insensitive to Cp, confirming that C3 
can be tuned independently of C^. (We ignore the vector meson because, as noted 
above, at a = 0.4 fm and P/V = 0.7 it is too heavy for its dispersion relation to 
be a reliable indicator of the Lorentz- violating errors.) 

From these results it is clear that for the D234c action, as for SW, the hadron 
dispersion relation is more sensitive to finite-a errors than the hadron mass. At 
P/V = 0.7 on an isotropic 0.4 fm lattice we get a perfectly sensible vector mass 
(Fig. 1^), but the dispersion relation is beginning to break down due to (9((am)^) 
errors. An anisotropic lattice would extend D234c's range to much higher quark 
masses at the same spatial lattice spacing, because the errors due to the quark 
mass are of order (atm)^. 



4.3 Redundant terms 

As a final test of the D234c action, we varied the parameter r. The r-dependent 
terms in the D234c action are, by design, approximately redundant, and therefore 
varying r should have little effect on the spectrum. Furthermore the extent to 
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Figure 5: Speed of light squared using the D234c action with couphng C3 = 1.2 
for pseudoscalar (P) and vector (V) mesons, at P/V — 0.76 (upper figure) and 
P/V = 0.6 (lower figure). 
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Action 


Wilson 


sw 


D234c 


a (fm) 


0.4 


0.4 


0.25 


0.4 


P 


16.0(7)% 


8(2)% 


1.0(1)% 


2.0(7)% 


V 


15.2(11)% 


8(2)% 


1.3(2)% 


2.0(5)% 


N 


9.0(14)% 


5(2)% 


0.7(2)% 


1.1(6)% 


D 


9.0(20)% 


8(2)% 


0.9(2)% 


1.5(6)% 



Table 4: Percentage change in hadron masses when the redundant coefficient r 
is changed from 1 to 2/3. Scaling violation due to redundant terms at r = 1 is 
therefore expected to be roughly 3 times this. 

which this is true gives us an indication of the size of the finite-a errors associated 
with the r terms, including the dominant correction, the clover term. Our results 
are summarized in Table ^. As expected, the variation with r is much smaller 
with D234c than with the SW or Wilson actions — the r terms are more highly 
corrected in the former. This "r-test" suggests that r-dependent finite-a errors in 
the static masses are of order 3-4% at 0.25 fm, and 6% at 0.4 fm for D234c. Our 
discussion above shows that these errors could be almost entirely due to radiative 
corrections to the clover coupling Cp- 

5 Conclusions 

Our studies show that at tree level the mean-link tadpole-improved D234c action 
is accurate to about 10% in hadronic masses, even for lattice spacings as large 
as 0.4 fm and meson masses as large as 2/a~l GeV. The meson dispersion relation 
(c^(p)) is similarly accurate out to three-momenta of order 1.5/a^750 MeV. It is 
markedly superior to the SW action in the limit of large masses or momenta. 

Our analysis indicates that the leading source of error in the hadron masses 
is the radiative correction, beyond tadpole-improved tree-level, to the clover 
coefficient Cp-, which enters in the 0{a) correction. Radiative corrections to other 
terms contribute probably less than 5% at a = 0.4 fm, while uncorrected errors 
are probably no more than a few percent. Thus the errors in the hadron spectrum 
in our 0.4 fm simulation can almost certainly be reduced to less than 10% by 
correcting Cp- 

We estimated the size of the quantum corrections to C^? by tuning Cp until our 
hadron masses, computed with a = 0.25 fm and 0.4 fm, agreed with a quadratic 
a = extrapolation of SCRI's (tadpole-improved) SW data. We found that 
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the quantum corrections to are of the right size to be perturbative. The 
perturbative prediction for Cp is not yet known, but we may compare with the 
clover coefficient in the SW action which, for Wilson's (unimproved) gluon action, 
has the perturbative expansion 1+0.46 as after mean-link tadpole-improvement [^, 
which agrees with non-perturbative estimates to within a reasonably sized 
error. 

We also demonstrated how to nonperturbatively tune the leading 
(9(a^) correction, the C3 term in the D234c action, by using restoration of 
rotational invariance. The correction to tadpole- improved tree-level is small, and 
is unimportant at the few percent level for a < 0.4 fm. 

As expected, all errors were much smaller at a = 0.25 fm, and radiative 
corrections were much less important. With a properly tuned clover coefficient, 
Cp, D234c should give (static and kinetic) hadron masses accurate to within a 
couple of percent at this relatively large lattice spacing. 

In the near future we will complete perturbative calculations of Ci;' and C3 
to compare with our nonperturbative tunings; nonperturbative tuning of Ci;' using 



PCAC restoration is another possibility |17|. We are also examining currents, form 



factors and other quantities with the D234c action. Perhaps most importantly, we 
are continuing our previous work on D234c for anisotropic lattices. This simulation 
technology is applicable to a range of phenomenological applications (hadronic 
spectra and decay constants, relativistic heavy-quark physics, high-momentum 
form factors), and we anticipate being able to obtain useful results with spatial 
lattice spacings in the range 0.2-0.4 fm. 
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Appendix A Galilean Quarkonium 

As explained in Section ^ Galilean quarkonium is a convenient fictitious state for 
the evaluation of ratios of lattice spacings. It is simulated using the conventional 
NRQCD quark action, without relativistic corrections: 

G(f,t + .) = (l - og)"(l - of )c/!,(l - «f ) (l - a^ycm (A.1) 
where 

a (A<'))^ a- 
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For details of the lattice derivatives and fields, see [0. 

For this paper, we needed the ratios of our two lattice spacings to each 
other, and to SCRI's (3 = 7.4 improved glue. We therefore performed Galilean 
charmonium simulations using all 3 types of glue, at a range of quark masses. The 
quark mass is characterized by the dimensionless quantity Mkin/(P — S). M^in is 
the kinetic mass of the quarkonium S-state, obtained from measurements of the 
energy of the S'-state at rest and with the smallest momentum on the lattice. This 
quantity is 6.8 for charmonium, so our Galilean quarkonium states were mostly 
lighter than charmonium. 

Our raw measurements of the P — S splitting are given in table ^ and the 
resultant lattice spacing ratios are in table |. We see that there is no dependence 
of the ratios on the quark mass within statistical errors, so that finite-a errors in 
our ratios are negligible. 

Using a potential model for Galilean quarkonium, we find that the radius of the 
meson changes by 25 — 30% over this range of quark masses, so any uncorrected 

errors would double in size. Interestingly, the potential model gives similar 
sensitivities of P — splitting to quark mass to those measured on the lattice (last 
line of table |^). 
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Mkin 









-s) 




p 


3.6 


4.0 


5.2 


7.2 


1.157 


1.0846(33) 


1.0545(60) 


1.0124(70) 


0.9667(80) 


1.719 


0.6863(35) 


0.6729(50) 


0.6317(80) 


0.5990(51) 


7.4 (SCRI) 


0.6638(36) 


0.6513(67) 






dln{P-S) 
dM^J{P - S) 


-0.06(1) 


-0.05(1) 


-0.04(2) 


-0.02(1) 



Table 5: Measurements of a{P — S) for Galilean quarkonium at various kinetic 
quarkonium masses Mkin. Also (final line), sensitivity of P — jS" to quark mass. 



a(/5= 1.157) 
a(/5= 1.719) 

a(/5= 1.157) 
a(SCRI,/5 = 7.4) 

a(/5= 1.719) 
a(SCRI,/5 = 7.4) 



jP-S) 

3.6 4.0 5.2 7.2 

0.633(4) 0.638(6) 0.624(9) 0.620(8) 
0.612(4) 0.618(7) 
0.967(7) 0.968(12) 



Table 6: Ratios of lattice spacings calculated using Galilean quarkonium, for a 
variety of quarkonium masses. Our glue {j3 = 1.157, 1.719) is mean-link tadpole- 
improved, SGRFs (/3 = 7.4) is plaquette tadpole-improved. 
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Appendix B Data Tables 





lattice 




a{P - S) 


a (fm) 


(MeV) 


1.157 


5^ X 10 


2.75 


0.925(7) 


0.400(4) 


495(4) 


1.719 


8^ 


1.80 


0.578(10) 


0.249(5) 


790(10) 



Tabic 7: Determination of lattice spacing from NRQCD calculation of spin- 
averaged charmonium P — S splitting, P — S = 458 MeV. The lattice size and 



bare quark mass 


. are given. 












D234c, am, 


,uark = 0.640 


SW, am. 


luark = 0.594 








73 


8=^ 


6=^ 


73 


8^ 


P 
V 
N 
D 


0.901(5) 
1.293(9) 
1.950(20) 
2.210(20) 


0.917(6) 
1.312(12) 

1.955(15) 
2.210(20) 


0.908(4) 
1.304(8) 

1.950(25) 
2.220(20) 


0.939(6) 
1.332(8) 
1.960(40) 
2.250(20) 


0.955(6) 
1.350(15) 

1.990(20) 
2.220(40) 


0.945(4) 
1.333(8) 

1.958(30) 
2.230(25) 



Table 8: Finite volume errors: Hadron masses in lattice units at P/V ^ 0.7, 
a = 0.25 fm, measured on 6^ x 20, 7^ x 20, and 8^ x 20 lattices. 
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a = 0.25 fm 


a = 0.40 fm 




D234c 


SW 


D234c 


SW 


quark 


0.68 


0.63 


1.002 


1.012 


P/V 


0.756(5) 


0.762(6) 


0.756(4) 


0.763(3) 


P 


1.051(4) 


1.085(4) 


1.583(3) 


1.497(7) 


V 


1.391(7) 


1.424(7) 


2.095(10) 


1.960(10) 


N 


2.100(30) 


2.150(24) 


3.220(15) 


3.000(20) 


D 


2.340(30) 


2.370(20) 


3.540(30) 


3.260(20) 



Table 9: Hadron masses in lattice units at P/V ~ 0.76. For a = 0.40 fm we used 
a 5^ X 18 lattice, for a = 0.25 fm an 8^ x 20. Physical spatial volume is the same. 





a = 0.25 fm 


a = 0.40 fm 




D234c 


SW 


D234c 


SW 


quark 


0.640 


0.594 


0.943 


0.970 


P/V 


0.694(3) 


0.708(8) 


0.703(4) 


0.699(8) 


P 


0.902(3) 


0.945(4) 


1.371(3) 


1.313(15) 


V 


1.300(4) 


1.333(8) 


1.950(10) 


1.878(10) 


N 


1.925(25) 


1.958(30) 


2.960(20) 


2.850(20) 


D 


2.170(20) 


2.230(25) 


3.380(20) 


3.150(30) 



Table 10: Hadron masses in lattice units at P/V ~ 0.70. For a = 0.40 fm we 
used a 5^ X 18 lattice, for a — 0.25 fm an 8^ x 20. Physical spatial volume is the 
same. 
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C [1^) 




c [^V ) 




mom 


ri 1 

Oi? — 1 


n 7o 

L^F — U. / Z 


r< 1 

Oi? — i 


O — U . ( z 












100 


0.94(4) 


0.91(3) 


1.08(3) 


0.91(3) 


101 


0.94(4) 


0.96(2) 


1.16(4) 


0.96(4) 


002 


0.90(2) 


0.93(4) 






no 


0.90(3) 


0.89(5) 


1.30(10) 


0.92(6) 


111 


0.84(5) 


0.98(9) 







Table 11: Sensitivity of to Cp: Hadron at P/F = 0.76, 5^ x 8 x 18, a = 0.4 fm 
lattice. The vector particle's (? shows definite sensitivity to a change in Cp {P/V 
kept constant) 











c\V) 






mom 




C?, = 1 


G = 1.2 




Cy, = 1 


C-A = 1.2 


001 


0.632(24) 


0.96(2) 


1.01(2) 


0.421(16) 


1.05(5) 


1.23(5) 


100 


0.550(13) 


0.93(3) 


1.01(2) 


0.342(13) 


1.08(3) 


1.23(6) 


101 


0.560(14) 


0.93(3) 


1.02(2) 


0.355(14) 


1.17(4) 


1.32(6) 


002 


0.510(16) 


0.90(2) 


0.99(2) 








no 


0.516(20) 


0.90(3) 


1.08(3) 


0.345(13) 


1.28(9) 




111 


0.505(60) 


0.84(5) 


1.01(8) 









Table 12: Sensitivity of to C3, the coeff of A^^^. Tree-level value is C3 = 1, SW 
has C3 = 0. The table shows hadron at P/V = 0.76 on a 5^ x 8 x 18, a = 0.4 fm 
lattice. Prom the pseudoscalar c^, C3 = 1.2 appears to be the results that would 
follow from non-perturbative tuning. 
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c^(P) 




c\V) 




mom 


SW 


D234c 


SW 


D234c 


100 


0.79(3) 


0.99(3) 


0.68(6) 


0.97(5) 


110 


0.76(4) 


0.95(4) 


0.65(3) 


0.96(4) 


111 


0.84(9) 


0.98(9) 






200 


0.70(9) 


0.94(8) 







Table 13: Hadron at P/V = 0.76, on an 8^ x 20, a = 0.25 fm lattice. 
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